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Abstract 

^ ■ By considering the continuum scaling limit of the A4 RSOS lattice model of Andrews- 

Baxter-Forrester with integrable boundaries, we derive excited state TBA equations 
describing the boundary flows of the tricritical Ising model. Fixing the bulk weights 
to their critical values, the integrable boundary weights admit a parameter £ which 
ON . plays the role of the perturbing boundary field 931.3 and induces the renormalization 

group flow between boundary fixed points. The boundary TBA equations determining 
the RG flows are derived in the 2) — ¥ ^(2,1) example. The induced map between 
distinct Virasoro characters of the theory are specified in terms of distribution of zeros 
of the double row transfer matrix. 



%*• 1 Introduction 

Quantum Field Theories with a nontrivial boundary have received recently a lot of atten- 
tion, due to their applications in Condensed Matter, Solid State Physics and String Theory 
(D-branes). A problem of great interest is the Renormalization Group flow between differ- 
ent boundary fixed points of a CFT that remains conformal in the bulk. Many interesting 
results have been achieved, and flows have been studied for minimal models of Virasoro 
algebra and for c = 1 CFT (see e.g. [T] and references therein). Numerical scaling functions 
for the flow of states interpolating two different boundary conditions can be systematically 
explored by use of the Truncated Conformal Space Approach [21 E]- 

A beta function can be defined for the boundary deformations, much the same as for 
the bulk perturbations of conformal field theories (Ij. The conformal boundary conditions 
can thus play the role of ultraviolet (UV) and infrared (IR) points of the flow. One gets 
out of the UV point by perturbing with a relevant boundary operator and gets into an IR 
one attracted by irrelevant boundary operators. 

Among the possible boundary perturbations of a CFT there are some that keep an 
infinite number of conservation laws. These are referred to as integrable boundary per- 
turbations. In this particular case experience suggest that there should be some exact 

1 Talk given by F. Ravanini at the 6th Landau Workshop "CFT and Integrable Models" - Chernogolovka 
(Moscow), Russia, 15-21 September 2002 
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method of investigation, based on transfer matrix and Bethe ansatz techniques. One of 
the most celebrated of these methods is the Thermodynamic Bethe Ansatz (TBA) IE] 
giving a set of non-linear coupled integral equations governing the scaling functions along 
the RG flow. In the following we illustrate how one can obtain TBA equations from a 
lattice construction for the integrable boundary flows. 

The simplest nontrivial example to consider is the second model in the Virasoro minimal 
series, the Tricritical Ising Model (TIM) with central charge c = j^. It has interesting 
applications in Solid State Physics and Statistical Physics. Its Kac table is given here 
below: 
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We have drawn only half of the table, by taking into account the well known Z 2 sym- 
metry (r, s) = (4 — r, 5 — s) that reduces the independent chiral primaries to 6 in the 
TIM. 

Conformal boundary conditions for minimal models with diagonal modular invariant 
partition function of type (A pj A p i) are classified to be in one to one correspondence with 
chiral primary fields in the Kac table. So there are 6 different conformal boundary condi- 
tions B( r>s) , r = 1, 2, 3 and s = 1, 2 for the TIM (A 3 , A 4 ). 

Let us consider in the two-dimensional plane (x, t) the TIM defined on a strip of thick- 
ness L, i.e in the region < x < L and — oo < t < +oo, with non-trivial b.c. -B( r ,s) 011 
the left (x = 0) and S( r / )S /) on the right (x = L) edge respectively. This situation will be 
denoted B( r , s )\(r',s,)- In the following we are interested in boundaries of type (r, s)|(l, 1). We 
shall denote for short B(r, s )\(i,i) = &(r,s)- For this particular class of boundary conditions 
the partition function simply reduces to one single character [Zj 

z B {r>s) {q) = X(r.s){q) 

where X(r,s)(<z) denotes the character of the irreducible Virasoro representation labeled by 
(r, s) at central charge c = j^. 
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For each given boundary condition #( r ,s) there is a set of boundary operators 4>( u ,v) that 
live on the edge. If we want to keep the b.c. unchanged all along the edge, these operators 
must be restricted to the conformal families appearing in the OPE fusion of the Virasoro 
family (r,s) with itself: (u,v) G (r,s) x (r,s). They are distinguished in terms of their 
conformal dimensions h r ^ s as relevant (h r , s < 1) and irrelevant (h r , s > 1). Of course only 
relevant perturbations break scale invariance at the boundary in such a way to get out of 
the fixed boundary point of a specific boundary condition and flow to another one. So 
if we want to consider possible relevant boundary perturbations of the TIM, i.e. QFT's 
described by the action 

/+oo 
dt(fi( u ,v)(x = 0,t) 
-oo 

- where S(r,s) denotes the action of TIM with boundary condition #( r , s ) - we have to restrict 
to the following possibilities 



boundary condition 


boundary perturbations 


#(1,1) = #(3,4) 




none 


#(2,1) = #(2,4) 





none 


#(3,1) = #(1,4) 


+ 


none 


#(1,2) = #(3,3) 


-0 


0(1,3) 


#(2,2) = #(2,3) 


d 


0(1,2), 0(1,3) 


#(3,2) = #(1,3) 


0+ 


0(1,3) 



Actually there are two physically different flows for each "pure" perturbation (i.e. con- 
taining only one operator 0( i2 ) or 0(i3)), flowing to two possibly different IR destinies. This 
is achieved by taking different signs in the coupling constant in the case of 0(i3) flows, and 
real or purely imaginary coupling constant in case of 0(i2) ones. The boundary condition 
#(2,2) can be perturbed by any linear combination of the fields 0(i,2) and 0(i,3). The symbols 
+, — , 0+, —0, 0, d represent other ways to denote the TIM conformal boundary conditions 
present in the literature 0121, and we report them here only to facilitate translations to 
the reader. 

Along the flow, going from UV to IR, the boundary entropy associated to each #( r , s ) 



9(r,s) = 




decreases Therefore we expect to have flows only between boundary conditions where 
the starting conformal boundary entropy is higher than the final one. 

The possible conformal boundary conditions have been studied extensively by Chim [8J 
and the flows connecting them by Affleck [H]. The picture can be summarized as in fig. 
HJ Integrability can be investigated in a manner similar to the bulk perturbations and it 
turns out that the flows generated by pure 0(i,3) and 0(i,2) perturbations are integrable. 
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Instead, the flow starting at -8(2,2) as a perturbation which is a linear combination of 0(i,2) 
and 0(1,3) is strongly suspected to be non-integrable. 

Notice the Z 2 symmetry of fig. [U which is strictly related to the supersymmetry of the 
TIM. Investigation of the supersymmetric aspects of the TIM boundary flows is, however, 
out of the scope of the present talk. 

We deduce TBA equations for the integrable flows from the functional relations for the 
transfer matrix of a lattice RSOS realization of the model. We will focus, as an example, 
on the flow obtained by perturbing the -8(1,2) boundary TIM by the relevant boundary 
operator ^(1,3) 

UV = B (1(2) ^ fi (2>1) = IR, Xl>2 (q) _> X2>1 (q) (1.1) 

This flow induces a map between Virasoro characters Xr,s(q) of the theory, where q is the 
modular parameter. The physical direction of the flow from the ultraviolet (UV) to the 
infrared (IR) is given by the relevant perturbations and is consistent with the g-theorem |3]. 
The approach outlined here, however, is quite general and should apply, for example, to 
all integrable boundary flows of minimal models. 

2 Scaling of Critical A4 Model with Boundary Fields 

It is well known that the TIM is obtained as the continuum scaling limit of the generalized 
hard square model of Baxter [TO] on the Regime III/IV critical line, i.e. the At RSOS 
integrable lattice model of Andrews-Baxter- Forrester [TT], with critical Boltzmann weights 

d c 



u 



a b 

/ d c \ sin(f-«) sinw /sinf sinf V /2 

\a b J sin | sinf Ism^sm^y 

satisfying the Yang-Baxter equation. The indices a,b,c,d = 1, ...4 take values on the 
A4 Dynkin diagrams, whose adjacency dictates which are the nonzero elements of the 
Boltzmann weights. 

In the presence of integrable boundaries, one must also introduce the so called boundary 
Boltzmann weights, satisfying the Boundary Yang-Baxter equations. Their general form 
has been given in [12], for our purposes here we select 
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All these Boltzmann weights are periodic under u — > w+7r. The scaling energies of the TIM 
are obtained jI3j from the scaling limit of the eigenvalues of the commuting double-row 
transfer matrices D(u) [H] 




of the A4 lattice model with N faces in a row. Integrability means [D(u),D(u)] = 0, 
Vu, v e C. The physical region is given by < u < |; crossing symmetry D(u) = D(| — w) 
implies that the isotropic point is u = j^. The matrix D(u) satisfies the property 

dim B (aL) = ( AN )^n 

where the matrix A is the incidence matrix of the A 4 Dynkin diagram. Here the notation 
D(a L) a R ) means the sub-matrix of D for fixed ol,<2r boundary spins. It is convenient to 
introduce the normalized double row transfer matrix 

t(u)=7>(u,&)D(u) 



6 



that satisfies the following functional equation: 



t(u)t (u + |) = 1 +t (u + . (2.1) 

For the interested reader, the expression for the normalization factor V(u, £l) is given in 
[T3] . Even if the equation for t looks boundary independent, the solutions are fixed by the 
behavior of the zeros, that are boundary dependent. The transfer matrix D(u) is entire 
and has only zeros; t(u) has instead two complex conjugate poles (coming from P(w,£l)) 
in the strip — ^ < Imu < We shall see that they will affect the boundary driving term 
in the TBA equations we are going to deduce. 

In a lattice of N horizontal and M vertical sites, with periodic boundary conditions on 
the vertical direction and cil,cir boundary conditions on the horizontal one, the partition 
function of the tricritical hard square model can be written as 

Z MN (u) = Ti n (D [aL , aR) (u) M ) ~ e- NMf ™M- Mf "»'MZ(q) 

where fbuik and f sur f are the bulk and surface contributions to the free energy respectively. 
Z(q) is the universal conformal partition function, with modular parameter 

Therefore the behavior of each transfer matrix eigenvalue is 

logd(u) = -Nf Mk (u) - f S urf(u) - ^ (A - ^ + nj sin 5m + o (J-^j 

with n a nonnegative integer. It is then logical to consider for the eigenvalues t(u) of t(u), 
the decomposition for N — > oo 

t(u) = f(u)g(u)l(u) 

where f(u) ~ O(N), g{u) ~ 0(1) and l{u) ~ 0(1/N). 

In general the conformal partition function expands over characters of the c = ^ 
Virasoro algebra 

Z{q) = Y^MtvXrM 

r,s 

where J\f r>s are nonnegative integers describing the multiplicities of Virasoro representations 
in the Hilbert space. 

Bauer and Saleur [11] used Coulomb gas techniques to show that if one chooses the 
boundary conditions leading to the configuration depicted in figEJ the single character 
X(r,s), i-e. the partition function of TIM with S(>,s) boundary conditions, can be realized 
in the scaling limit N — ► oo. The configuration can be built by taking the boundary 
Boltzmann weights described above, with = r, = s and Re^,L = jq- This latter 
choice kills the weight alternating r and r — 1 on the left at the isotropic point u = f^. 
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Figure 2: Lattice boundary configuration leading to Cardy partition function 



At criticality there exist integrable boundary conditions associated with each conformal 
boundary condition (r, s) carrying 0,1 or 2 arbitrary boundary parameters £l,£r [T2j . To 
consider boundary flows we must fix the bulk at its critical point and vary the boundary 
parameters £l,£r. However, these "fields" are irrelevant in the sense that they do not change 
the scaling energies if they are real and restricted to appropriate intervals jTH]. Explicitly, 
the cylinder partition functions obtained [EH E21 from (r, s) integrable boundaries with 
£,l,£,r real are independent of and given by single Virasoro characters 



Z{q) = XrM- 



(2.2) 



The reason for this is that in the lattice model the fields £l and £r control the location of 
zeros of eigenvalues of D(u) on the real axis in the complex plane of the spectral parameter 
u whereas only the zeros in the scaling regime, a distance zlogiV out from the real axis, 
contribute in the scaling limit N — > oo. The solution is to scale the imaginary part of 
£l,£r as logiV. This is allowed because £l,£,r are arbitrary complex fields. 

We find that only one parameter at the time (we call it £l) is needed to induce the 
^(1,3) flows and is identified as the thermal boundary field: 



£l ~ ^(1,3)- 

For the flow (|1.1|) . following [1 3j , we consider the (r, s) = (2, 
boundary on the right (with no boundary field) and the (2, 1) 
boundary field We scale £l as 



(2.3) 

1) boundary with the (1, 1) 
boundary on the left with 



Re(6 



-7T. 



Im(& 



-e + iogiv 



(2.4) 



10 ' ^ ' 5 

where £ is real. In terms of boundary weights, the (1,2) boundary is reproduced in the 
limit Im(£z,) — > ±oo whereas the (2, 1) boundary is reproduced in the limit Im(£i) — > 0. 
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The scaling behavior of the At lattice model at the boundary fixed points is described 
by the TBA equations of O'Brien, Pearce and Warnaar (TH] (OPW). We generalize their 
analysis to include the scaling boundary field ()2.4JI . 



3 Classification of zeros 

The zeros of D(a L ,a R )( u ) can be classified in two strips, the first given by — ^ < Re(u) < || 
and the second by ^ < Re(tt) < -jr. We refer to them as strip I and II respectively. Each 
zero has to be accompanied by its complex conjugate, and there are only two zeros on 
the real axis depending on Re£z, that do not play any role in the following. So we are 
led to consider only zeros in the upper half u complex plane excluding the real axis. For 
N sufficiently large general considerations supported by numerical inspection allow to say 
that there exists 

• mi 1-strings at ^5 + iy^ and n\ 2-strings at — i^?_±__LlZL _j_ with j = 0, 1 in strip I 

• m 2 1-strings at ^ + iy^' and n 2 2-strings at ( 2j + 1 ) 7r + with j = 0, 1 in strip II 

The precise classification can be done in terms of the relative ordering (with respect to 
their imaginary parts) of the 1 and 2-strings. This can be achieved by giving the sequence 
of numbers Ijy of 2-strings with > yf 1 . In the case of B^i >2 ) this must be supplemented 
[IH] by a parity of strings a — ±1. Because we order the zeros as y^. x > y%\ the quantum 
numbers must satisfy > I® > ... > I® > 0. Moreover, in the second strip the largest 
quantum number must be less or equal the number of 2-strings of type 2, n 2 > I®. The 
2-strings in strip I are infinite in number at the scaling limit (see below), and their quantum 
numbers are unbounded. We write this sequence as 

j«i/ (2) j (2) i^) with o-l ±l ioiB w 

1 — K 1 ! 5 J 2 5 •••5 i mi l J l ) 1 2 ) •••) 1 m 2 )^ wlul 1 for B{2 i) 



There are constraints on the numbers of strings, analyzed in detail in [T3], depending on 
the specific B(r,s) condition 

• for B(i t 2) ( m i odd, m 2 even) 

N + m 2 + a mi - a . 
n 1 = mi , n 2 = m 2 (3.1) 



for £>(2,i) {mi odd, m 2 odd) 

N + m 2 
rii = mi 



n 2 



mi + 1 



- m 2 



(3.2) 
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In the scaling limit N — >• oo the number of 2-strings in strip I becomes infinite: n\ ~ O(N), 
while the number of the other three types of roots (1-strings in both strips and 2-strings 
in strip II) remains of order 0(1). In [TH] it was shown that the pattern of zeros fixes the 
energies (here m = (7711,7712) and C is the A4 Cartan matrix) 

, 1 mCm mi+m 2 AA/:) 

1=1 fc=i 

t=l k=l 

Out of these energy expressions it is possible to reconstruct the characters xi,2 and \2,i- 
First consider the so called finitized characters [T6] 

M 

e v^o («)/ 

where is the largest eigenvalue of the double row transfer matrix with (r, s) boundary 
conditions and N sites. The limit N — > 00 of such objects can be shown to give the well 
known formulae for the TIM characters (for details see [T7]). 

Once the characters have been related in such a strict way to the distribution of zeros 
in strips I and II, it is also possible to see how they are mapped one another along the 
flow 13(i,2) ~^ 13(2,i). For this purpose it is more convenient to consider the reverse flow 
13(2,i) ^(1,2) from IR to UV. The counting of the total number of zeros in the upper u 
plane for a 6(2,1) configuration is N, while for 6(1,2) is N — 1. Therefore one zero must 
disappear to infinity along this reverse flow (or appear from infinity in the physical flow). 
We have found 3 mechanisms that seem to exhaust the possible ways for such a phenomenon 
to be realized [TJj 

• A. The top 1-string in strip 2 flows to +00, decoupling from the system while m 1? ni 
and 712 remain unchanged. This mechanism applies in the IR when £) = Iml = 
and produces frozen states in the UV with a — 1 and £) = 0. 

• B. The top 2-string in strip 1 and the top 1-string in strip 2 flow to +00 and a 
2-string comes in from +00 in strip 2 becoming the top 2-string. Consequently, each 
Ij 1 ^ decreases by 1 and each 1^ increases by 1. This mechanism applies in the IR 

when Iml > and Iml — and produces states in the UV with a = — 1 and either 
m 2 = or Iml > 0. 

• C. The top 2-string in strip 2 flows to +00 and a 1-string in strip 2 comes in from 

(2) 

+00. Consequently, each If, decreases by 1. This mechanism applies in the IR when 

(2) (2) 
Imi > and produces states in the UV with a = — 1 and I m2 = 0. 
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Observe that 

• these mappings are in fact one-to-one 

• the 1-strings in strip 1 are never involved 

• the change in parity of m 2 is consistent with (|3.1j) and (I3.2B . 

• for A,B: m% v = m} 2 R — 1 and for C: m% v = mt^" + 1. 

By doing such substitutions m^ v — > ml R ± 1 in the finitized characters, it is possible to 
show explicitly that Xipio) ~^ X2,i{q) under the physical flow #(1,2) — > £>(2,i)- See [TJj for 
details. 



4 Derivation of Boundary TBA Equations 
4.1 On the lattice 

The recurrence relation (|2.1J) for the O(N) piece yields the following functional equation 
for the function f(u) 

f(u)f(u +£) = 1 

Plugging its solution (for details see [THj) back into the functional equation (12.11) it yields an 
equation for g(u) that in turn can be solved in the first analyticity strip — ^ < Reu < y^. 
In the following, for any h{u) = f(u),g(u),l(u) we introduce the functions 

hi{x) = h ^— + i- j , h 2 {x) =h\— + i- 

each one being analytic in the corresponding strip I or II. 
The recurrence relation for g(u) yields 

2 x x + 5Im£ L x - 5Im£ L 

Qi(x, El) — tanh — tanh tanh 

yu -sl; 2 2 2 

and an analogous formula for g 2 (x) that we do not need in the following. The lattice TBA 
equations contain the boundary contribution encoded in the function gx(x). It becomes 
exactly 1 at the conformal points and takes £_l dependent values during the flow. Notice 
that gi(x) depends on ^ only through its imaginary part. 

Inserting this value for g\{x) into the recurrence equation again, one can determine 
a constraint to be satisfied by the last pieces U(x), in the form of an integral equation. 
Putting all these things together we arrive finally to two coupled integral equations to be 
satisfied by ti(x) and t 2 (x) 

mi 

-logti(x) = - log gi (x) + lo g tanh ^(.x-^p)- AT*log(l + t 2 ) (4- 1 ) 

3=1 

m2 ^ 

-\ogt 2 {x) = -A r log(tanh 2 -) + ^logtanh-(x-4 2) )-A'*log(l+ti) (4.2) 

2 2 

fc=i 



ii 



where the kernel of the convolution is: 



K(x) 



27rcosh(x) 

Observe that in there is an implicit N dependence that will take a role performing 
the scaling limit. Only the 1-strings appear explicitly in the TBA equations, v k being 
their imaginary parts. Their location is dictated by the quantization contitions 



-log (-tiO 
-log [~t 2 {v 



(2)< 
k , 



(1)' 
k , 



(2) • 



(1) • 

n k 7Ti 



where n k = 2(mi — k + I k L> ) + 1 — m 2 and 



0-h 



n 



(2) 



2(m 2 -A; + jf ) ) + l-m 1 . 



(4.3) 



(4.4) 



4.2 The scaling limit 

In the scaling limit N — > oo where Im^ rescales as in (|2.4j) we define hi(x) = liniTv^oo hi(x+ 
log A/"). We have 

&Or,0 =tanh^ (4.5) 

with £ — * — oo and £ — > +oo corresponding to S(i,2) and S(2,i) respectively. 
Eqs. (|4.3|4.4|) fix the following large N behavior for v®: 

y = lim (5v — log N) = finite val. 

N — > OO 



Performing this scaling limit on (|4.1l4.2J4.3J4.4j) . and introducing the so called pseudoen- 
ergies 

Ei{x) = - log (-U(x)) 

and the functions 

L, : (x) = \o K (l - e- £ * (x) ) 



we obtain the following final TBA equations: 

— log tanh — tarm I J g 



— K * L 2 (x) 



m 2 



5 2 (x) = 4e x — ^ log tanh 



(2) 

yl> -x 



K*LAx) 



k=l 



The definition of the scaling theory must be completed by the energy formula 

E(E) - — = - lim RE(R) = 



1 

7T 



£ 2 e-»S"-r^l„ g (l 

i= l •'-oo 



-£2(a;)> 



(4.6) 



(4.7) 
(4.8) 



(4.9) 
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where the limit R — * means that all this computations are performed at the critical 
temperature of the TIM. The boundary parameter £ appears only implicitly in (|4.9|) . The 
last piece required to complete the scaling picture is the location of the zeros. From l|4.3J4.4J) 
we obtain: 

£2(2/!^ — 2—) = n^Vi first strip zeros (4-10) 
J 2 

(2) ^" f2) 

£ i(^fe 2 = * second strip zeros (4-11) 

The two-string zeros satisfy similar equations, using the substitution y^' — > 2^ . 
An immediate consequence is that Li(x) vanish on the zeros of the same strip: 

U{ V f) = 0. (4.12) 

The inverse statement is not true, in the sense that there can be points where L, vanishes 
without corresponding to some zero. 

We do not deal here with the numerical analysis of the TBA equations found above. 
Details can be found in |T2j and will be an important part of our paper in preparation [T^l 
where we explore the other integrable flows of the TIM. Generalizations to higher minimal 
models and to other rational CFT's should follow the same lines, although they can become 
more technically involved. 
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